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ABSTRACT 

We provide new, high-resolution A{K S ) extinction maps of the heavily reddened Galactic midplane 
based on the Rayleigh- Jeans Color Excess ("RJCE") method. RJCE determines star- by-star red- 
dening based on a combination of near- and mid-infrared photometry. The new RJCE-generated 
maps have 2' x 2! pixels and span some of the most severely extinguished regions of the Galaxy — 
those covered with Spitzer+IRAC imaging by the GLIMPSE-I, -II, -3D, and Vela-Carina surveys, 
from 256° < / < 65° and, in general, for \b\ < 1°-1.5° (extending up to \b\ < 4° in the bulge). Us- 
ing RJCE extinction measurements, we generate dereddened color-magnitude diagrams and, in turn, 
create maps based on main sequence, red clump, and red giant star tracers, each probing different 
distances and thereby providing coarse three-dimensional information on the relative placement of 
dust cloud structures. The maps generated from red giant stars, which reach to ~18-20 kpc, probe 
beyond most of the Milky Way extinction in most directions and provide close to a "total Galactic 
extinction" map — at minimum they provide high angular resolution maps of lower limits on A(K S ). 
Because these maps are generated directly from measurements of reddening by the very dust being 
mapped, rather than inferred on the basis of some less direct means, they are likely the most accurate 
to date for charting in detail the highly patchy differential extinction in the Galactic midplane. We 
provide downloadable FITS files and an IDL tool for retrieving extinction values for any line of sight 
within our mapped regions. 

Subject headings: Galaxy: structure - Galaxy: disk - reddening, dust, infrared techniques 



1. INTRODUCTION 

The existence of interstellar dust and its light- 
obscuring effects, though postulated more than a cen- 
tury and a half ago by Struve (1847) as one (incor- 
rect) solution to Olber's Paradox, was not definitively 
proven until Trumpler's (1930) insightful experiment to 
compare distances of open clusters derived in separate 
extinction-dependent and extinction-independent meth- 
ods. Though the amount of extinction resulting from this 
dust was historically often approximated with a stan- 
dard accumulation by distance under an assumed plane- 
parallel distribution (i.e., the cosecant law; e.g., Shane 
& Wirtanen 1967, Peterson 1970, Noonan 1971, Sandage 
1973, de Vaucouleurs, de Vaucouleurs & Corwin 1976, 
de Vaucouleurs & Buta 1983, Fabbri et al. 1986; but see 
Kron & Guetter 1973, Holmberg 1974, Knude 1996 for 
examples of discussions of simple variations about the 
basic cosecant law), it was always well-understood that 
the distribution of Galactic dust is in fact quite irregu- 
lar]]] Once it is appreciated that an irregular selective 
obscuration of a more regular distribution of stars, and 
not an irregular distribution of the stars themselves, is 
shaping the net configuration of observed starlight, then 
any observation of the density of — or light from — stars 
in the Galactic plane of the Milky Way clearly illustrates 
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1 Even relatively recently, Knude (1996) went through great 
pains to determine whether the cosecant law can at least provide 
a reasonable approximation in the North Galactic polar cap and 
concluded that only a piecewise application of the cosecant law 
with varying coefficients can work, and only in certain ranges of 
Galactic latitude and longitude. 



this patchiness (e.g., Herschel's [1785] drawings of the 
distribution of stars in the Milky Way, or the early pho- 
tographs of the Galactic plane by Barnard et al. 1927). 

For extragalactic studies, a long-standing default strat- 
egy has been simply to steer clear of the "Zone of Avoid- 
ance" at low Galactic latitudes altogether, rather than 
contend with assessing the amount of foreground extinc- 
tion. Unfortunately, Nature is not always so accommo- 
dating with the placement of important and interesting 
sources, and corrections for foreground dust — gener- 
ally attributable to dust in our own Milky Way — are 
almost always needed; even at the highest Galactic lat- 
itudes, one can encounter appreciable amounts of dust 
(e.g., the North Galactic Pole has A(B) as high as 0.05 
to 0.1 mags, depending on the reddening calibration and 
extinction law used; e.g., Hilditch et al. 1983, Knude 
1996, Peek & Graves 2010). 

Early attempts to map the variation of total line-of- 
sight Galactic extinction explored a number of tech- 
niques. For example, under an assumed constancy of 
dust-to-gas ratio, the net column density of HI, easily 
measured via 21-cm observations, could be converted into 
a dust column density (e.g., Lilley 1955, Sturch 1969, 
Savage & Jenkins 1972, Knapp & Kerr 1974, Burstein & 
Heiles 1982). Under the assumption of a homogeneous 
spatial distribution of galaxies, the counts of galaxies to 
a given magnitude limit were also used to infer the degree 
of extinction along a line of sight (Burstein & Heiles 1978, 
1982). Similarly, the assumption of spatial constancy of 
the stellar distribution (at least on relatively small angu- 
lar scales) has been employed to derive extinction using 
star counts — comparing stellar densities across a clearly 
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extinguished region (often on the edges of a dark cloud) 
to those in nearby "control" regions (e.g., Wolf 1923, Bok 
1956, Froebrich et al. 2005). 

Directly measuring the effects of dust via such affected 
tracers is the most reliable method to map the Galactic 
extinction. Extracting the reddening and extinction ef- 
fects of foreground dust from the intrinsic spectral energy 
distributions (SEDs) of observed sources has tradition- 
ally required direct, detailed spectroscopic information 
of each source to classify it explicitly and thereby infer 
its intrinsic properties independent of observed colors — 
quite a demanding prospect in most situations. In recent 
decades, a large body of work has amassed that attempts 
to determine the total amount of extinction foreground 
to heavily-reddened sources by photometric means; un- 
der the assumption of a homogeneous color distribution 
of either stars or galaxies (or particular classes of either) , 
deviations from the mean observed colors are attributed 
to interstellar reddening. The decreased temperature de- 
pendencies of infrared (IR) colors, combined with the de- 
creased extinction effects at IR wavelengths, has led to 
the wide adoption of IR and optical-IR color excesses as 
high-quality gauges of interstellar reddening, especially 
in highly extinguished regions (since e.g., Jones et al. 
1980, 1984, Smith 1987). Lada et al. (1994) introduced 
a powerful and efficient method of extinction mapping 
with the NICE (Near Infrared Color Excess) technique, 
which applies elements of both IR excess measurements 
and star-counting techniques to large sets of IR photom- 
etry; this method has undergone significant evolution 
(NICER, Lombardi & Alves 2001; V-NICE, Gosling et 
al. 2009; NICEST, Lombardi 2009) and has been ap- 
plied to a series of dense clouds and dark nebulae (e.g., 
Alves et al. 1998, Lada et al. 2009, Lombardi et al. 2010). 
The color excess method has also been applied to galaxy 
colors (e.g., Peek & Graves 2010), though galaxy ensem- 
bles are obviously best-suited to gauging total line-of- 
sight Galactic reddening. Greater accuracy may be at- 
tained through use of objects with well-defined (and lim- 
ited ranges of) colors like RR Lyrae stars, globular clus- 
ters, or red clump stars. Another method advocated by 
Burstein and collaborators (Burstein et al. 1988, Faber et 
al. 1989) is based on the observation that the integrated 
spectra of the stellar populations of giant elliptical galax- 
ies are strikingly uniform, so that the intrinsic (B — V)q 
colors of these systems can be inferred on the basis of the 
spectroscopic measurement of their Mg2 line index. 

Early investigations into the consistency of extinction 
derived with these various methods generally led to dis- 
appointing results, with no simple relationships found 
amongst the different dust proxies (Heiles 1976). Partly 
this is because of the failure of the very assumptions on 
which each method is based — e.g., the high "noise" in 
galaxy counts due to large scale structure, or variations 
in the dust-to-gas ratio (Seki 1973, Jenkins & Savage 
1974, Bohlin 1975, Savage et al. 1977, Burstein & Heiles 
1978, Fong et al. 1987, Burstein 2003). 

Currently the most commonly used large-scale red- 
dening maps are those by Schlegel, Finkbeiner & Davis 
(1998, "SFD" hereafter), which have been derived from 
100/im emission maps using the COBE/DIRBE and 
IRAS/ISSA data. These maps have proven reliable for 
differential line-of-sight reddening in low extinction re- 
gions, albeit with some uncertainties in zero-points (see 



discussion by Burstein 2003). Unfortunately, these maps 
have a resolution of only ~6', which is particularly prob- 
lematical where the reddening varies dramatically on 
much smaller angular scales. In addition, a number of 
studies have shown that the SFD maps overestimate 
the extinction at low latitudes by as much as 50% in 
A(V) (Stanek 1998a,b, Arce & Goodman 1999a, Chen 
et al. 1999, Ivans et al. 1999, von Braun & Mateo 2001, 
Choloniewski & Valentijn 2003, Rocha-Pinto et al. 2004, 
Dutra et al. 2003a,b, Cambresy et al. 2005, Amores & 
Lepine 2005, 2007, Yasuda et al. 2007, Peek & Graves 
2010). A physical reason for the complications of the 
SFD maps at low latitudes is explored in the analysis in 
Paper I — namely that far-infrared- emitting dust, which 
provides the basis of the SFD extinction maps, is sim- 
ply not a good proxy for the dust causing extinction 
at shorter wavelengths. Comparison to other extinction 
tracers reveals that many of the SFD "high extinction" 
knots actually trace supernova remnant or star-forming 
environments, and a number of cold, dense clouds can 
be identified in extinction that lack the far-IR emission 
actually traced by the SFD maps. 

Thus, the midplane extinction in the Milky Way is far 
from being unambiguously mapped, and for studies of 
the Milky Way itself, the patchy intermingling of dust 
and stars remains a particularly troublesome obstacle. 
Both photometric and spectroscopic parallax determina- 
tions must include an accurate assessment of the extinc- 
tion foreground to a star to gauge its distance. While 
studies of specific dark clouds using stellar IR color ex- 
cess methods (particularly the extensive NICE family of 
techniques) have achieved remarkable success at mapping 
extinction, reliable assessments of extinction foreground 
to stars and astronomical sources not associated with 
these dark clouds have proven more elusive. 

Now, more accurate measurements of individual stellar 
extinction levels, including full three-dimensional spatial 
distribution information, is becoming possible through 
combined near- and mid-IR observations of the Rayleigh- 
Jeans flux distributions of stars (Majewski et al. 2011, 
Paper I hereafter); because the long wavelength SEDs 
of normal stars have approximately the same shape 
(i.e., colors) on their Rayleigh- Jeans tails, the observed 
Rayleigh- Jeans colors contain information on the fore- 
ground reddening to a star explicitly. We show in Paper 
I that this Rayleigh- Jeans Color Excess (RJCE) method 
is superior to the SFD method for mapping the extinc- 
tion by midplane dust and is an improvement over the 
other IR color-excess methods due to the improved ho- 
mogeneity of RJCE's adopted color indices. An addi- 
tional advantage of the RJCE method over its predeces- 
sors is that the former cleanly separates reddening from 
intrinsic stellar colors and preserves the latter for use 
in the identification of stellar types (e.g., main sequence 
versus red clump versus red giant); thus, we can use the 
RJCE method to fairly reliably isolate the most distant 
stellar proxies (i.e., red giants) to create the best assess- 
ment of the total reddening along a given line of sight. 
Alternatively, one can isolate a rather reliable standard 
candle (in the form of red clump stars) to bring greater 
accuracy to distance estimates. In combination, stellar 
types of different intrinsic luminosity can yield the dis- 
tance distribution of the dust along the line of sight. 

In this contribution to the present series of papers on 
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the RJCE method, we make new, two-dimensional ex- 
tinction maps of the Galactic midplane across the entire 
Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire (GLIMPSE-I; Benjamin et al. 2003; GLIMPSE- 
II/-3D Churchwell et al. 2009), and Vela-Carina (Spitzer 
PID 40791) IRAC survey areas. Section 2] summarizes 
how the maps were created; the reader is directed to Pa- 
per I for a fuller description and analysis of the RJCE 
method and its application to different color index com- 
binations and different stellar types. In Section [3] we 
present a full atlas of these maps based on use of main 
sequence, red clump, and red giant tracers, and in Sec- 
tion [4] we alert the reader to various caveats for using the 
maps. Electronic versions of these maps have been made 
publicly available, and we describe how to access them 
in Section [U 

2. CREATION OF THE MAPS 
2.1. Summary of the RJCE Method 

IR imaging, like that in 2MASS (Skrutskie et al. 2006) 
and UKIDSS (Lucas et al. 2008), can more efficiently 
penetrate the obscuration in the dusty midplane and 
Galactic Center than observing at shorter wavelengths, 
and much work has been done using IR observations of 
low-latitude stars as probes of Galactic structure and the 
interstellar medium (ISM) there. Measuring "color ex- 
cesses" is one common use of this NIR photometry — 
comparing the observed colors of a star (or the average 
of a group of stars) to the expected intrinsic colors to de- 
termine the amount of line-of-sight reddening, which can 
then be converted into an extinction under an assumed 
extinction law. Figure 1 in Paper I demonstrates some of 
the primary complications when this method is applied 
to broadband NIR colors in isolation. First, the spread 
in intrinsic NIR color is significant: A, F, G, and K stars 
span about a magnitude in ( J — H)o and (J — K s )o and 
almost half of a magnitude in (H — K s )q. In addition, the 
NIR reddening vector very closely parallels the intrinsic 
stellar locus, leading to high degeneracy in derived red- 
dening estimates based on observed stellar colors. 

Both of these critical problems are greatly amelio- 
rated by using longer- wavelength, longer-baseline colors 
including MIR photometry from, e.g., the Spitzer Space 
Telescope's IRAC instrument. These longer- wavelength 
filters sample the Rayleigh- Jeans tail of a normal star's 
SED; the slope of this tail (and thus the colors measured 
from it) has very little dependence on stellar temper- 
ature, so the intrinsic variation in MIR and NIR-MIR 
colors is significantly smaller than in NIR colors. For 
example, the A-K stars mentioned earlier differ by only 
~0.1 mag in (H — [4.5/i])o- Furthermore, the reddening 
vector in these colors is very nearly perpendicular to the 
stellar RGB locu^J facilitating more accurate estimates 
of the stellar reddening. The RJCE method takes ad- 
vantage of these NIR-MIR color properties to produce 
reliable estimates of the line-of-sight reddening to each 
individual star. Analysis of many factors, including in- 
trinsic color variation and typical photometric quality 
(see Paper I), has led us to adopt (H — [4.5//]) as the 
best single NIR-MIR color for stellar dereddening. 

2 It is true that the reddening vector more closely parallels the 
locus of cool main-sequence stars; however, these are very sparse 
in the 2MASS and GLIMPSE surveys 



Also shown in Paper I are comparisons of RJCE- 
derived extinction maps to those of SFD as well as to 
maps of 13 CO emission, 8/xm dust and PAH emission, 
and NIR star counts. The RJCE maps bear a striking 
structural similarity to the 13 CO and star count maps, 
both expected to be good tracers of cold, dense extinc- 
tion regions, while the SFD maps more closely resemble 
the [8.0//] images. 

Moreover, the RJCE method directly measures the red- 
dening on a star-by-star basis, and is therefore superior 
to methods that use proxies to calculate extinction (e.g., 
dust emission, HI, or galaxy counts). The RJCE method 
relies only on the availability of stars, which are, fortu- 
nately, very abundant in those parts of the Galaxy where 
the effects of dust are strongest. 

2.2. Defining the Stellar Samples 

Arguably, the RJCE method's most powerful and 
unique feature is the ability to recover a star's approxi- 
mate stellar type after the reddening has been estimated. 
The strong temperature dependence of the corrected NIR 
colors (especially those involving J-band) allows us to 
derive an approximate photometric temperature (hence 
spectral type and, based on the apparent magnitude, 
the luminosity class as welQ. For the initial publicly- 
available maps presented here, we have used this powerful 
feature to split the 2MASS+ Spitzer stellar sample into 
three subsamples: main sequence (MS) dwarfs, which are 
dominant for 0.0 < ( J — K s )o < 0.5, red clump (RC) gi- 
ants for 0.55 < (J — K s )o < 0.85, and red giant branch 
(RGB) giants for 0.85 < (J-K 3 ) < 1.4, where (J-K 3 ) 
for each star is ca lcula ted using the RJCE method de- 
scribed in Section |2.1j and Paper I. See Figure [I] for an 
example of a corrected color-magnitude diagram(CMD) 
and the rationale for the adopted stellar color divisions. 
Each of these stellar subsam ples probes a different range 
of heliocentric distance (see £4.2 and Figure [5]), so maps 
using each tracer, taken together, provide coarse three- 
dimensional information about the midplane dust distri- 
bution, based solely on which of the actual stellar tracers 
are being reddened. 

3. EXTINCTION MAPS OF THE MILKY WAY MIDPLANE 

Figure [3] contains the "Version 1" RJCE extinction 
maps, in units of A(K S ), over ~170° of longitude in the 
Galactic midplane. For \l\ < 65°, the latitude range is 
centered on the Galactic equator (and includes \b\ < 1° 
for most longitudes but up to \b\ < 4° in the bulge). In 
the longitude range covered by the Vela-Carina survey 
(295° < I < 255°), the map center varies between b = 
+0.5° and b = —0.5°, but the latitude coverage is con- 
sistently 2° tall. 

The basic procedure used to make the maps is as fol- 
lows: Using the reliable band-merged 2MASS+IRAC 
Catalog] we selected stars having detections and pho- 
tometric errors <0.2 mag in all four of the J, H, K s and 

3 There is certainly dwarf/giant degeneracy at redder NIR colors, 
but Galactic stellar population models (e.g., TRILEGAL; Girardi 
et al. 2005) guide us as to which magnitudes, for a given color, 
we may expect a low fraction of late- type dwarf contaminants. 
Fortunately, these models show that late-type dwarfs are in any 
case very rare in the 2MASS+GLIMPSE catalog. 

4 |http : //irsa . ipac . caltec h . edu/data/SPITZER/GLIMPSE/doc/ 
glimpsel_dataprod_v2 . . pdf 




Fig. 1. — Three color-magnitude diagrams (CMDs) of a 4 deg 2 field at (l,b) = (42.0°,0°) for all sources with magnitude errors less 
than 0.5 mag in the three 2MASS and four Spitzer-IRAC bands. Sources are binned in bins of 0.025 in J — K s and 0.04 in K s . (a) The 
raw 2MASS data. (6) The RJCE-corrected CMD. (c) An extinction-free model of the same field from TRILEGAL (Girardi et al. 2005) 
with photometric errors; color-coded to highlight the different stellar populations (MS-main sequence, RC-red clump, and RGB-red giant 
branch). The vertical dashed lines in (6) and (c) indicate the color divisions between MS, RC and RGB stars. 

[4.5/i] bands. The extinction was calculated for each star 
separately, using 



A(K 3 ) = 0.918 x(H- [4.5/i] - 0.08) 



(1) 



from Paper I. Each star was dereddened and classified 
as MS, RC, or RGB using the RJCE method and the 
stellar color bins described in Section [2j Each of these 
three subsamples was binned separately into a spatial 
map with 2' resolution, and two extinction maps were 
created for each subsample (for a total of six distinct 
maps) using the median and 90th percentile values of 
A{K S ) measured within each spatial pixel. Similar maps 
were created for the larger sample of all stars, with no 
color cuts. All of these extinction maps were Gaussian 
smoothed (FWHM=1.4 / ) to reduce noise. Not only do 
the 90th percentile maps represent the majority of the 
extinction affecting each stellar subsample, they are also 
the least likely to be biased by contaminants from, e.g., 
foreground cool dwarfs (in the case of the RGB maps) 
or photometric errors. All extinction values are given in 
magnitudes of A(K S ). We also provide star count and 
A{K S ) dispersion (not to be confused with the uncer- 
tainty) maps for each subsample. Pixels with less than 
two stars in them are set to white in the A{K S ) quick- 
look images (black in the starcount quick-look images) 
and have a value of 999999.0 in the queriable FITS ta- 
bles. 

Finally, we also include "lower limit total extinction" 
maps produced using only the IRAC [3.6/i] and [4.5/i] 
bands; these are 90th percentile extinction maps made 
using all stars with detections and photomet ric e rrors 
<0.5 mag in both of these bands (see Section 4.1). To 
derive extinctions from ([3.6/i]- [4.5/i]) we use, 

A(K S ) = 9.188 x ([3.6/i] - [4.5/i] + 0.235), (2) 

which was found by linear fitting the E(H—[4.5/j])- 
derived A{Ks) values and the ([3.6/i] — [4.5/i]) colors for 
all stars in the longitude range 20° < I < 30° and pho- 
tometric errors in H, [3.6/i], and [4.5/i] less than 0.07 
mag (see Fig. |2|. The intrinsic scatter and photomet- 
ric errors in (pTo/i] — [4.5/i]) are generally higher than for 



(H— [4.5/i]), and, therefore, the derived extinctions are 
generally less reliable and mainly useful for high extinc- 
tion regions (where the [H— 4.5/i] maps often drop out). 
Including the ([3.6/i] — [4.5/i]) extinction maps, there are 
a total of 18 maps in the downloadable FITS files, as 
listed in Table [I] Examples of the extinction maps are 
also shown in Figure [3] 
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Fig. 2. — Linear fitting of E(H — [4.5/i]) -derived A(K S ) values as 
a function of ([3.6/i] — [4.5/i]) color for all stars with photometric 
errors in the H, [3.6/i], [4.5/i] bands less than 0.07 mag. The values 
were binned in bins of 0.25 mag in A(K S ) (red dots) and the robust 
linear-fit perfomed with the binned values. The best-fit line is 
shown in black. 
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Fig. 3. — A subset of the publicly available maps. The first three maps are 90th percentile A(K S ) extinction maps as derived with the 
RJCE method using (H — [4.5//]) colors and (a) MSTO, (b) RC, and (c) RGB stars. Panel (d) shows the 90th percentile A(K S ) extinction 
map derived using [3.6/x]-[4.5/x] colors of all stars (note the different colorscale ranges compared to panels a-c). Finally, panel (e) shows the 
starcounts map of all stars. The GLIMPSE-I data cover 10.0 < I < 65.4° and 295.0 < I < 350.3° and -1 < b < +1°, while the Vela-Carina 
data cover 255.3 < I < 295.0° and —1.5 < b < +0.5°. Pixels with no stars are shown as white in the extinction maps and black in the 
starcounts map. 
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TABLE 1 
Maps in Downloadable FITS Files 



HDU a 


Pop. Type 


Map Type b 


Statistic 


1 


MSTO 


RJCE (H - 


'4.5/i' 


) Extinction 


Median 


2 


MSTO 


RJCE (H - 


'4.5/x' 


) Extinction 


90th percentile 
Dispersion 


3 


MSTO 


RJCE (H - 


4.5/i 


) Extinction 


4 


MSTO 


Starcounts 




5 


RC 


RJCE (H - 


'4.5/x 


) Extinction 


Median 


6 


RC 


RJCE (H - 


4.5^ 


) Extinction 


90th percentile 


7 


RC 


RJCE (H - 


4.5/i 


) Extinction 


Dispersion 


8 


RC 


Starcounts 




9 


RGB 


RJCE (H - 


4.5/i 


) Extinction 


Median 


10 


RGB 


RJCE (H - 


4.5/i 


) Extinction 


90th percentile 


11 


RGB 


RJCE (H - 


4.5/i 


) Extinction 


Dispersion 


12 


RGB 


Starcounts 




13 


ALL 


RJCE (H - 


4.5/i 


) Extinction 


Median 


14 


ALL 


RJCE (H - 


4.5/i 


) Extinction 


90th percentile 


15 


ALL 


RJCE (H - 


4.5/i 


) Extinction 


Dispersion 


16 


ALL 


Starcounts 




17 


ALL 


RJCE [3.6/i]- 


-[4.5/i] Extinction 


Median 


18 


ALL 


RJCE [3.6/i]- 


-[4.5/i] Extinction 


90th percentile 



1 Extension in the FITS file. 

2 All extinction maps are in magnitudes of A(K S ). 

3 Not to be confused with an uncertainty estimate. 
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We note that due to the lower density of stars in the 
outer Galaxy data, there are a larger number of "empty" 
pixels in the RGB extinction maps in these directions 
than in the more populated inner Galaxy (although the 
RC and MS extinction maps suffer much less from this 
deficiency). Therefore, the RGB maps at these longi- 
tudes should be used with caution. Fortunately, our es- 
timates of th e distance range covered by the RC maps 
(Section |4~2| ) suggest that in these sparser outer Galaxy 
regions, the farthest RC stars are very likely to probe ex- 
tinction to the edge of the Galactic disk; thus, the 90th 
percentile A(K S ) RC maps may be used in place of the 
RGB maps where the latter are unhelpfully incomplete. 

3.1. Three Dimensional Information 

The derived extinction maps provide coarse three- 
dimensional information on the intrinsic interstellar dust 
distribution. Because each of the three stellar subsam- 



ples probes different distance ranges (see Section 4.2), 



comparison of the maps for a given section of sky revea 
information about the distance, as well as the thickness, 
of individual clouds. The /=50-60° panels of Figure 3] 
demonstrate this capability. For example, the ~0.5° wide 
cloud centered near (Z, b) ~ (54. 8°, 0.8°) in the RGB map 
(Figure |3fc) appears in neither the MS nor the RC map, 
which indicates that it is a distant feature, lying be yond 
the reach of the RC stellar sample (~8 kpc; Section 4.2) 



but in front of the RGB sample. Similarly, the filamen- 
tary structure near (Z, b) ~ (59.0°, 0.5°) is visible in the 
RC and RGB maps but not the MS one, which suggests it 
is closer than the former cloud. An even closer structure 
lies at (Z, b) ~ (53. 7°, 0.5°), as evidenced by the fact that 
it appears in the relatively short-range MS map as well 
as the more distant ones. In future work, we will be able 
to refine these cloud distance measurements after more 
careful calculation of the stellar distance distribution. 

3.2. Uncertainty Estimates 

Factors contributing to the uncertainty of RJCE ex- 
tinction measurements are discussed extensively in Pa- 
per I. In summary, the uncertainty sources of (1) intrin- 
sic scatter in the stellar colors, (2) observed photometric 
uncertainties, and (3) uncertainty and variation in the 
extinction law combine to yield a typical individual stel- 
lar extinction uncertainty of <0.11 mag in K s . The use of 
median and 90th percentile statistics in the maps further 
reduce t he extinction uncertainty p er pix el, as described 
i n, e.g ., Lombardi & Alves ( |2001| ) and Dobashi et al. 
([2009). 



4. CAVEATS TO USAGE OF THE MAPS 
4.1. Limited Photometry in Highly Reddened Regions 

Apparent in many of the extinction maps in Figure [3] 
are sharp "holes," generally at the higher-latitude edges 
of outer disk fields or at the center of high-extinction 
regions. The former are mostly due to an intrinsic lower 
stellar density (i.e., a true lack of stars) and could be 
remedied by reducing the map resolution; currently, for 
simplicity, we mask these pixels as white (see Section [3]). 

The holes at high extinction, however, are artifacts 
due to photometric limits in the adopted catalogs (here, 
the 2MASS PSC and GLIMPSE/ Vela-Carina; see Paper 
I for discussion of the merged catalog limits), where a 



particularly dense cloud has a total extinction sufficient 
to increase the magnitude of background stars beyond 
the limit of the catalog. When this cloud is relatively 
nearby, the result is an almost total lack of stars with 
which to measure the extinction; when the cloud is far- 
ther away, the measured extinction in the affected pixels 
is dominated by less-reddened foreground stars, so that 
the A(K S ) value of the central extinction "peak" is re- 
ported to be smaller than in the rest of the cloud. 

Automatically identifying which "holes" are due to 
which cause (s) is a nontrivial task. For the purposes of 
mapping extinction, there are three approaches to deal- 
ing with the second problem (i.e., artificially low A(K S ) 
values due to the presence of a nearby dense cloud): (1) 
acquire deeper photometry, (2) ignore (identify and flag) 
the data in such pixels, or (3) use photometric bands with 
even lower dust sensitivity in an attempt to "see" beyond 
the cloud. Option (1) will become increasingly feasible 
as large ongoing and future IR surveys release their data, 
but until then, deeper IR photometry is limited to much 
smaller patches of sky (mostly targeted at dark clouds or 
the central Galactic regions), not amenable to large-scale 
or diffuse extinction mapping. We have chosen option (2) 
for the extinction maps presented here, but we enable op- 
tion (3) by including ([3.6/i] — [4.5/x]) extinction maps for 
possible use in high extinction regions. 

Our analysis of catalog magnitude limits in Paper I 
suggests that the use of MIR data alone (i.e., not requir- 
ing a JHK S detection) holds great promise for detecting 
stars behind extremely dense clouds. The extinction at 
3.6/im and 4.5/im is only ~0.04xt4(V) (Indebetouw et 
al. 2005, Zasowski et al. 2009), compared to -0.2 x A(V) 
for the H band (Cardelli et al. 1989). Figure E] demon- 
strates that by using IRAC's [3.6/i] and [4.5/1] bands 
alone one can trace not only the extinction peaks appear- 
ing as "holes" in the primary RJCE (H— [4.5/i]) map but 
also an overall higher level of extinction. This is because 
these redder MIR bands are less susceptible to dust red- 
dening than the NIR bands and so probe farther, while 
detections in JHK S , the lack of which is the dominant 
cause for highly-reddened stars to be extinguished out 
of the band-merged 2MASS+ Spitzer catalog, are not re- 
quired. However, as discussed in greater depth in Paper 
I, a number of factors — including the greater intrin- 
sic scatter in ([3.6/i] — [4.5/i] )o, the increased photomet- 
ric uncertainties at longer wavelengths, and the inability 
to distinguish stellar type (and hence identify the most 
distant extinction probes) because all stars have essen- 
tially the same ([3.6/i] -[4.5/i]) color — render IRAC-only 
colors less suited for producing dereddened NIR CMDs 
or three-dimensional dust maps. Nevertheless, in our 
publicly- available maps (Section [3|, we include those de- 
rived using E{ [3.6/i] — [4.5/i]) as an option for users in- 
terested in estimates of the total line-of-sight extinction. 

4.2. Distance Limits 

Obviously, t he c atalog magnitude limits discussed 
above (Section |4.1| ) translate to distance limits in the 
maps and define the typical depth to which each map 
probes. Characterization of these limits is crucial to un- 
derstanding how to interpret the extinction data con- 
tained in the maps. This characterization was done in 
two ways. 

First, we analyzed the drop in stellar density at the 



8 



Nidever, Zasowski, & Majewski 




Fig. 3. (cont.) 



RJCE Extinction Maps 




(cont.) 



10 



Nidever, Zasowski, & Majewski 
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Fig. 4. — Total extinction maps (in units of A{K S ) mags) made by the RJCE method, but with use of the (a) (H — [4.5/x]) colors and 
(b) ([3.6/x] — [4.5/i]) colors as the selective extinction measure (where we have adopted ([3.6/x] — [4.5/x])o = —0.1 as a baseline color). In 
both cases, the stars used are those satisfying a 0.5 mag photometric uncertainty limit in only the two filters used in the color-excess 
measurement. Panel (c) shows the difference between the maps in panels (a) and (b). The few white pixels indicate bins containing <2 
stars. 

faint edge of the RGB in dereddened CMDs throughout 
the GLIMPSE and Vela-Carina survey coverage. In sum- 
mary, we compared the faint edge of the RGB to a suite 
of stellar isochrones spanning — 1.3 < [Fe/H] < +1.7 dex 
and 9.0 < log(Age) < 10.05 yrs (Girardi et al. 2002), 
leaving distance as a free parameter. This provided us 
with a rough estimate of the distance of the farthest stars 
in each CMD, as well as an estimate of how likely those 
stars were to be at the magnitude limit of the catalog (as 
opposed to representing a genuine drop in stellar density, 
such as at the edge of the Galactic disk) . The procedure 
used here is very similar to the one described in Paper I 
with the main difference that the disk edge (at a given 
color bin) was determined by finding the magnitude for 
which 90 percent of all stars in that bin are brighter. 
This updated method was found to produce reliable fits 
to the faint edge of the RGB and the results are pre- 
sented in Figure [5] The filled red circles indicate the 
derived Galactic positions for the most distant stars in 
the Al = 2° wide fields that were fit reasonably well with 
isochrones. Figure [5] and our comparison of the CMDs to 
the catalog magnitude limits are consistent in suggesting 
that our extinction maps from the most distant tracers 
(i.e., those made using the mean and 90th percentile of 



RGB extinction; see Section |3| likely measure the true 
total Galactic extinction. 

• MSTO: We selected all TRILEGAL model stars in 
a test field at (Z,6)=(30°, 0°) in the CMD region 
0.0 < (J - K s ) < 0.4 and 12.5 < K Sj0 < 13.5, 
where a "clump" of MSTO stars is seen at these 
Galactic coordinates. The median distance of these 
simulated stars is ~2 kpc and the 90th percentile 
is ~4 kpc. Therefore, the MSTO distances are 
roughly 2-4 kpc. 

• RC: Because red clump stars are standard candles 
(almost independent of [Fe/H]; Salaris & Girardi 
2002) and their absolute magnitude is known fairly 
accurately (Mk s = —1.54 ± 0.04; Groenewegen 
2008), we were able to calculate distances directly 
for our RC sample at all longitudes. Stars in the 
RC color range (0.55 < ( J — K s )o < 0.85) were se- 
lected and binned in 3' pixels. Then, for each pixel, 
the RC star with the 90th percentile extinction was 
selected and its distance was calculated using its 
dereddened K s magnitude and Mk s = —1.54. Fi- 
nally, the median distance was calculated in longi- 
tude bins of 2°; these are shown as a solid, jagged, 
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Fig. 5. — Distances of the stellar tracer populations used for the extinction maps as a function of Galactic longitude. Red filled circles 
i ndica te the maximum distance probed by the red giant branch as determined from isochrone fits to the lower edge of the dereddened RGB 
(j|4.2) . Median RC distances (in 2° bins) are shown by the black line, while the typical main sequence distance of ~3 kpc is shown in blue. 
1 he solid circles indicate Galactocentric radii of 12 and 15 kpc. The shaded region indicates the longitude range covered by our extinction 
maps. 

black line in Figure [5) These distances are roughly above are on average 15 kpc. Therefore, the RGB 

6-10 kpc with an average distance of 8.4 kpc. distances are roughly 15-24 kpc. 



• RGB: Both a TRILEGAL model and isochrones 
were used to estimate the RGB distances in a 
test field at (Z,6)=(30°, 0°). Sources with 0.85 < 
(J — K s )o < 1.4 and K Sj q < 13.5 were selected 
from the TRILEGAL model as RGB stars (only 
four red dwarfs contaminated this sample). The 
median distance for these sources is ~15 kpc and 
the 90th percentile is ~19 kpc. In addition, a 5 
Gyr, solar metallicity Padova isochrone was used 
to obtain distances for observed stars with 29° < 
I < 31°, -1° < b < +1° and RGB colors of 
0.85 < (J — K s )o < 1.25 (since the isochrone ter- 
minates at 1.25). Distances were derived assum- 
ing these were RGB stars. The median isochrone 
distance is ~10 kpc, the 90th percentile distance is 
~24 kpc, and the median value in a 90th percentile 
binned distance map (with 3' pixels) is 22.9 kpc. In 
comparison, the "CMD edge" distances described 



• ALL: Since the vast majority of the stars in the 
data set are red clump stars, the distance limits 
for maps that make use of all stars are likely to be 
similar to the limits for RC stars. 

The reader is referred to Figure [5j which shows dis- 
tances of the stellar tracers. The RGB maps probe to 
the edge of the Galactic disk and can be considered close 
to "total line-of-sight extinction" maps in nearly all di- 
rections covered by the maps. 

5. PUBLIC ACCESS TO THE MAPS 

To encourage usage, additional comparisons, and feed- 
back from the astronomical community, we have made 
public 515 deg 2 of midplane and bulge extinction maps 
using the RJCE method outlined in Paper I and the 
present paper — i.e., using the (H — [4.5//]) color with 
iJ-band photometry from 2MASS and [4.5/i] data from 
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Spitzer-IRAC . In addition, we provide E([3.6jj] — [4.5//])- 
derived extinction maps for "all stars" which are use- 
ful in heavily-extinctecl regions, and 2MASS star count 
maps that are roughly anti-correlated with extinction. 
The regions included in this public release are those 
covered by the GLIMPSE-I, -II, -3D, and Vela-Carina 
surveys, spanning 255.3 < I < 65.4° and, in general, 
for \b\ < 1° — 1.5° (the Vela-Carina data is centered 
slightly off the midplane to track the highest extinc- 
tion regions) but up to b ~ 4° for certain longitudes 
(for GLIMPSE-II and -3D). These maps are available at 
http://www.astro.virginia.edu/rjce/ as quick-look 
GIF images and downloadable FITS files along with 
querying software written in IDL. 

Included on the website are instructions for the down- 
load and use of the maps and software, as well as notes 
regarding the known caveats. We welcome any feedback 
from the community that may help improve these tools, 
and we plan to post refined and improved maps as our 
methodology evolves. 

6. SUMMARY AND FUTURE WORK 

We have developed a method for measuring line-of- 
sight extinction that combines near- and mid-infrared 
photometry (from 2MASS + Spitzer-IRAC) to calculate 
stellar reddening in a process that is (i) robust to varia- 
tions in intrinsic stellar type and (ii) recovers the intrinsic 
stellar type. This Rayleigh- Jeans Color Excess (RJCE) 
technique represents a major improvement on previous 
color-excess measurements because it includes reddening 
corrections of temperature-sensitive near-infrared colors 
which allows us to create dereddened CMDs (extremely 
useful for stellar population studies) and select reason- 
ably pure samples of stars of specific stellar type or lu- 
minosity class. In Paper I of this series, we described the 
RJCE method in great detail. This installment presents 
a large, publicly- available atlas of two-dimensional, high 
angular resolution extinction maps created from RJCE- 
derived extinction values. 

We take advantage of RJCE's ability to discriminate 
stellar spectral type to produce extinction maps gen- 
erated independently from three groups of stars domi- 
nated by distinct stellar populations (main sequence [MS] 
dwarfs, red clump [RC] giants, and red giant branch 
[RGB] stars), each of which probes different heliocen- 
tric distance ranges and thus traces the net dust column 
to different distances. Comparison of these maps for any 
given section of the Galactic midplane provides rough 
three- dimensional information on the intrinsic dust dis- 
tribution and allows clouds and cloud complexes to be 
isolated to specific distance ranges. We present a to- 
tal of ten different types of extinction maps (in units 
of A[if s ]), each encompassing ~515 deg 2 of the dusti- 
est regions of the Galactic midplane, using data from 
the 2MASS-PSC merged with Spitzer-IRAC photometry 
from the GLIMPSE-I/-II/-3D and Vela-Carina surveys. 
For each of the three stellar tracer subsamples (MS, RC, 
and RGB) and "all stars", we provide maps giving two 
different statistical estimates of the (H — [4. 5 jj])- derived 
extinction: the median and 90th percentile A(K S ) in 
each spatial pixel. The three 90th percentile maps prob- 
ably provide the most reliable estimate of the total inte- 
grated extinction affecting each subsample population to 
the maximum distance probed. Additionally, we provide 



median and 90th percentile ([3.6//] — [4. 5/i])- derived ex- 
tinction maps of "all stars" , which provide reasonable ex- 
tinction estimates for very dusty regions where the other 
maps fail due to lack of stars. 

Section [4] provides caveats to the applicability of our 
maps as well as distance limits for the different tracer 
populations. We find that our RGB 90th percentile maps 
very likely represent the total Galactic extinction in near 
all directions. Our 90th percentile RC maps probe to 
an average distance of ~8 kpc, while the MS maps are 
limited to 2-4 kpc. 

The possibility of using MIR-only colors to estimate 
extinction where NIR-MIR colors are unavailable (due to 
the increased extinction in the NIR bands or the lower 
resolution of 2MASS) is also explored. Though MIR-only 
colors do not carry enough information to allow stellar 
type selection (and thus reliable distance estimates), we 
find that they can provide extinction estimates in the 
most highly extinguished portions of our maps (such as 
cloud cores and regions near the Galactic center). Thus, 
we include in our public atlas extinction maps derived 
from MIR-only reddening to serve as strong lower limits 
to the total extinction; this is most useful in the inner 
Galaxy, where even the RGB maps from typical RJCE 
colors are insufficient to probe to the far edge of the 
Galactic disk. 

An extensive program for application and improvement 
of the RJCE methodology is anticipated. We are working 
to improve the RJCE extinction estimates using a more 
refined treatment of stellar properties, and the ongoing 
APOGEE survey (Majewski et al. 2010) will allow us to 
compare RJCE-derived stellar types to spectroscopically- 
determined ones. We are evaluating the benefit of using 
optimized "hybrid" -color schemes (i.e., merging extinc- 
tion estimated using long-baseline NIR-MIR colors and 
MIR-only colors) where high extinction or confusion lim- 
its the depth of the NIR photometry. Our ultimate goal 
is to apply this dereddening methodology to ascertain 
the three-dimensional distribution of both dust and stars 
within the Galactic midplane. We have already taken 
advantage of our extensive data set to characterize the 
wavelength dependence of IR extinction throughout the 
midplane (Zasowski et al. 2009), and work is under way 
to use the new, reliably-cleaned CMDs to explore heavily- 
reddened or otherwise elusive Galactic structures, such 
as the central bar(s) and the outer spiral arms and the 
edges of the disk. 
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